function y = computeFaceCenter(f)
nodes = f.getNodes();
l = nodes.length;
cx = 0;
cy = 0;
for i=1:l
    X = nodes(i).getPosition();
    cx = cx + X(1);
    cy = cy + X(2);
end
y = [cx/l cy/l];
